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Abstract 

We compute the cross section for exclusive production of a photon associated with a heavy 
quarkonium H of charge-conjugation parity C = +1 from e+e^ annihilation into a virtual photon at 
the center-of-momentum energy ^/s = 10.58 GeV. The nonrelativistic QCD factorization formulas 
for the differential and total cross sections are obtained at leading order in the strong coupling 
and in the relative velocity of the heavy quark in the quarkonium rest frame. The predicted 
cross sections for the S'-wave spin-singlet cases are about 80 fb and 50 fb for H = rjc and ric{2S), 
respectively. Among P-wave spin-triplet charmonia, Xci has a particularly large cross section of 
about 14 fb. The cross sections for bottomonium states rjh and Xbj are about 3 fb. A rough estimate 
of the background reveals that the signal significances for charmonium processes are sufficiently 
large enough to be detected with ease with the integrated luminosities available at the present B 
factories. 

PACS numbers: 12.38.-t, 12.39. St, 13.66.Bc, 14.40. Gx 
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I. INTRODUCTION 



The measurements of the cross section for exclusive J/ip + rjc production from e'^e" 
annihilation carried out by the Belle [1, 2] and BABAR [3] collaborations have triggered 
rapid progress in the heavy-quarkonium theory based on the nonrelativistic QCD (NRQCD) 
factorization approach [4]. The original empirical cross section in Ref. [1] was greater than 
the theoretical predictions by Braaten and Lee [5] and by Liu, He, and Chao [6] by an order 
of magnitude. Scenarios that the signal may contain other final states like J/ip + J/ip [7, 8] 
or J/'^+glueball [9] were ruled out by an updated Belle analysis [2]. There was an argument 
that an appropriate choice of the light-cone wave function [10] may enhance the theoretical 
prediction significantly, which was disfavored by a following study [11]. Zhang, Gao, and 
Chao reported large corrections of next-to- leading order in the strong coupling ag [12] and 
the result was confirmed by Gong and Wang [13], but the corrections were not large enough 
to resolve the discrepancy between experiment and theory. In the mean time, Bodwin, 
Kang, and Lee introduced a new method of resumming relativistic corrections of a class 
of color-singlet contributions to all orders in v [14], where v is the relative velocity of the 
heavy quark or antiquark in the meson rest frame. The method was applied to compute 
the relativistic corrections to e'^e~ J/tp + rjc [15]. Finally, Bodwin, Lcc, and Yu reported 
that the discrepancy had been resolved within errors [16], after including both the resummed 
relativistic corrections and next-to-leading order QCD corrections by Zhang, Gao, and Chao. 
In addition, the first proof of the factorization theorem for exclusive charmonium production 
in B decay and exclusive two-charmonium production from e'^e~ annihilation to all orders 
in perturbation theory in QCD was presented very recently by Bodwin, Garcia i Tormo, and 
Lee [17]. 

As shown in Ref. [5], at leading order in v, a non-negligible amount of the cross section for 
the exclusive production of charge-conjugation parity C = +1 charmonium (H) associated 
with a J/iJj from e~^e~ annihilation into a virtual photon comes from the interference between 
QCD and QED subprocesses, where the QED subprocess involves 7* — > -|- 7* followed by 
the photon fragmentation 7* — > J/t/j. The interference contributions are about -|-30%, +5%, 
—6%, and 10% for H = rjc, Xco, Xa, and Xc2, respectively [5]. This means that the pure QED 
contributions are about 10%, 0.3%, 0.4%, and 1% of the total cross sections, respectively. 
We notice that, in comparison with the QED subprocess for e"'"e~ H + J/ip, the rate for 
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e+e~ — > if + 7 is enhanced by a scaling factor of 1/(0:1'^) ~ 10^ for a charmonium, where 
a. is the QED couphng. Based on this rough estimate, we expect that the cross section for 
e^e~ ^ if + 7 can be greater than that for e^e~ H + J /ip hy 2 orders of magnitude for 
H = rjc- For the P-wave spin-triplet charmonia, the enhancement factors range from about 
1 to 10. Therefore, as long as the severe background from e^e~ — > X + 7 is well controlled 
in the recoil mass {mx) regions near the H resonances, the processes e^e~ H + j can be 
detected by analyzing the photon energy spectrum in e'^e" — > X + 7. 

In this paper, we compute the cross section for the exclusive process e'^e" — > if + 7 at 
the center-of-momentum (CM) energy ^/s = 10.58 GeV, where if is a heavy quarkonium 
of charge-conjugation parity C = +1. We consider S'-wavc spin-singlct states rjc, ric{2S), 
and 7]b and P-wave spin-triplet states Xcj and Xbj, where J = 0, 1, and 2} Perturbative 
calculation is carried out at order a^a^ and at leading order in v. The remainder of this 
paper is organized as follows. In Sec. II, we compute the short-distance coefficients for the 
NRQCD factorization formulas for the e'^e~ — > if -|- 7 cross sections. The numerical values 
for the input parameters such as long-distance NRQCD matrix elements are given in Sec. III. 
The cross sections for e~^e~ ^ if -|- 7 at the B factories are predicted in Sec. IV, where a 
rough estimate on the background is also given. Then, we summarize our results in Sec. V. 

II. PERTURBATIVE CALCULATION 

In this section, we derive the NRQCD factorization formula for the exclusive process 
e+e" —> if -|- 7 at leading order in both and v, where if is a heavy quarkonium of charge- 
conjugation parity C = +1. The quarkonia if that we consider in this paper are the S'-wavc 
spin-singlet CSq) state, rig, and the P-wave spin-triplet (^Pj) states, xqJj for <^ = 0, 1, and 
2, where the heavy quark Q is the charm or the bottom quark. Because the C parity of 
the final state is —1, at leading order in a, an exclusive H + j final state can be produced 
from e'^e~ annihilation into a virtual photon. A heavy quarkonium with C — —1 can be 
produced associated with a photon from e'^e~ annihilation into two photons, which we do 
not consider here. 

^ Throughout this paper we suppress the identifier (IP) for any IP quarkonium. 
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FIG. 1: Feynman diagrams for the exclusive production of a C = +1 heavy quarkonium associated 
with a photon from e'^e~ annihilation into a virtual photon at leading order in a and a^. 

A. Amplitude 

The Feynman diagrams for the exclusive process e~^{k2) e'{ki) H{P, \h) + ■y{k, A) at 
order a^a^ are shown in Fig. 1. Here, ki, k2, P, and k are the momenta for the e~, e"^, H, 
and 7, respectively. The hehcities of the H and 7 are denoted by \h and A, respectively. 
The ^'-matrix element for the process is given by 



where — e is the electric charge of the electron and is the polarization four-vector of the 
photon. The leptonic current in Eq. (1) is defined by 



The factor A'^e*^, in Eq. (1) corresponds to the amplitude for 7*(Q) H{P, \h) + 'y{k, A), 
where Q is the momentum of the virtual photon 7*. 

B. Effective vertex 

A'^ in Eq. (1) contains nonperturbative contributions of the QQ pair that evolves into 
the heavy quarkonium H. In the NRQCD factorization approach [4], the hadronic tensor 
A'^ can be expanded as a linear combination of NRQCD matrix elements that involve the 
long-distance nature of the QQ pair. A corresponding short-distance coefficient of each 
NRQCD matrix element is insensitive to the long-distance nature of the pair and calculable 
perturbatively. According to the velocity scaling rules of NRQCD, NRQCD matrix elements 
are classified in powers of v [4]. Thus, the expansion can be truncated at a given order of 



iMh{Xh,X) = -2-L^^^^(AH)e;,(A), 



(1) 



= v{k2)'^^u{ki). 



(2) 
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accuracy. The velocity scaling factor of an NRQCD matrix element is determined by the 
spectroscopic state of the QQ pair, [QQ]i(^^'^^Lj), and that of the H into which the QQ 
pair evolves. Here, s, L, and J are the spin, orbital angular momentum, and total angular 
momentum of the QQ pair, respectively, and the subscript 1 indicates that the pair is in a 
color-singlet state. 

For the exclusive process 7* H + j, the color-singlet QQ pair with the spectro- 
scopic state identical to that of the H contributes at leading order in v. In order to de- 
termine the corresponding short-distance coefficient, we compute the QQ analog Aq" of the 
hadronic tensor A'^, where ^Q^(^*+^Lj)e*j, is the amplitude for the perturbative process 

7*(Q) [Qip) Q{p)]i(^''^~^^Lj) + j{k). Here, the four-momenta p and p for the Q and the 
Q, respectively, can be expressed in terms of the total momentum P and half their relative 
momentum q: 

P=^P + q, (3a) 
P = - (3b) 

where P • g = 0. In the QQ rest frame, P = {2E, 0), g = (0, qr), £; = (m^ + q2)V2^ and m is 
the mass of the heavy quark Q. At leading order inv, P can be identified as the quarkonium 
momentum. The amplitude for the perturbative process 7*(Q) — > Q{p)Q{p)+l{k, A) is given 
by 

Mni*{Q) ^ Q{p)Q{p) + i{K A)] = u{p)A^''v{j))e*^^{\). (4) 
The tensor A*^"" in Eq. (4) is defined by 

^z^- = e|e' [-i^Ki-p - k)Y + YHp + kh"] 1, (5) 

where eg is the fractional electric charge of the heavy quark Q, A(p) — {"^ -\- m) / {p^ — m^), 
and 1 is the unit matrix of the color SU(A^c) for Nc — 3. 

A^'' in Eq. (4) may contain various [Q{p)Q{p)\ii^^'^^Lj) contributions. The QQ analog 
A'^i^^'^^Lj) of the hadronic tensor can be obtained from A'^^ in Eq. (5) by restricting the 
QQ pair to have an appropriate spectroscopic state. A given spin state of the color-singlet 
QQ pair can be projected out from Eq. (4) by replacing the outer product v{p)u{p) with 
the corresponding spin projector given in Ref. [18]. The color-singlet spin-singlet projector 
is denoted by Hi and the color-singlet spin-triplet projector is by H^e*(As), where e(As) is 



5 



the polarization for the spin-triplet state with the helicity A^. Expressions for the Hi and 
rig are given by [18] 

^^^tf-mhy+2EW + n,)^ 1 ^^^^ 



4:^/2E{E + m) y/K' 



-.1. (6b) 



^ 4:V2E{E + m) ^ ^/Wc 

Wc still need to project out an appropriate orbital-angular-momentum state. In order to 
project out the 5"- wave contribution at leading order in one must put g = 0. By making 
use of Eqs. (5) and (6a), we obtain the QQ analog Aq(^Sq) of the hadronic tensor for the 
S-wsive color-singlet spin-singlet state: 

Ai^^CSo) = Tr[^'^'^ni] . (7) 

5=0 

where the trace is over both spin and color indices. 

The P-wave contribution at leading order in v is proportional to the first derivative of 
Eq. (4) with respect to q. By making use of Eqs. (5) and (6b), we obtain the QQ analog 
e*p{Xs)el{Xe)AQ''^ of the hadronic tensor for the P-wave color-singlet spin-triplet state, where 
e{Xe) is the polarization for the P-wave orbital-angular-momentum state and 



(8) 

9=0 



Like Eq. (7) the trace in Eq. (8) is over both spin and color indices. The QQ analog 
Aq'i^Pj, Xh) of the hadronic tensor for the [QQ]i{^Pj) for J = 0, 1, and 2 can be read off 
by projecting out the diagonal, antisymmetric, and symmetric traceless components with 
respect to the vector indices for s and L [5]: 

A'qCPo) = (9a) 
A^^^P,, Xh) = A^Q^'^^^e^rapP^'e'SiXH), (9b) 



2EV2 

A'Q''eP2,XH)=A^Q^'-e},^{XH), (9c) 



where Xh is the helicity of H. in Eq. (9b) and e^^ in Eq. (9c) are the spin-1 polarization 
vector for the [QQ]i{^Pi) pair and the spin-2 polarization tensor for the [QQ]i{^P2) pair, 
respectively. Note that = and en'^a — 0- Because we consider the color-singlet 
contributions at leading order in v, the polarization vector and tensor are treated to be 
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identical to those for Hi^Pi) and H(^P2), respectively. The tensor I^^, appearing in Eq. (9) 
is defined by 



P P 



{2EY' 

As a result, we obtain the effective vertex for the perturbative process 7' 

[gg]i('^+'i^j)+7as 



(10) 



Q^-P^ 

8V2ele\Q^ - 3P^) 
\fP^{Q'^ - P'^Y 

P2(g2 _ p2)2 



(11a) 
(lib) 



P^PaQpe'""^" - P^'PaQpe^'"" - \{Q^ - P')Pae'''""' 



A'^i'P, 



+^(g'-i^VV 



(lie) 



(lid) 



In the limit 0, the expressions in Eq. (11) reduce into the effective vertices for the 

decay [QQ]i{^'^^^Lj) — >■ 77. We observe that the vertex for ^Pi in Eq. (11c) is proportional 
to Q^. Because of this reason, 7* [QQ]i{^Pi) + 7 is not forbidden while the amplitude 
for [QQ\\{^Pi) — > 77 vanishes. 



C. Perturbative matching 

The QQ analogs Aqi^'^^^Lj) (11) for the hadronic tensors contain perturbative NRQCD 
matrix elements. At the squared amplitude level, those perturbative NRQCD matrix ele- 
ments at the leading order in v are (0|(!?[*^*^^^^ '^"'■'(^*'+^Lj)|0) for the color-singlet '^^^^Lj 
state, where Oi is the local four-quark operator of NRQCD [4]. In general, the production 
matrix elements are related to the NRQCD matrix elements for the decay under vacuum- 
saturation approximation, which is valid up to corrections of relative order [4]. For 
exclusive production through the electromagnetic interaction, the v-leading perturbative 
NRQCD matrix elements for the ^5*0 and ^Pj states are simply expressed at the amplitude 
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level as 



([Qg]i(^^o)|^^x|0) = 2^V2iVe, (12a) 

^([QQ]i(=^Po)|V'^ (-iB • a) x\0) = 2E^/2N:\q\, (12b) 

73([W]i('A)|^^ (-fB X o- ■ e^,) x|0) = 2£;^2iVe|q|, (12c) 

J^iiQQU'P^)]^^ (-iBV)e^l) x|0) = 2Ey/2N:\q\, (12d) 

where -0^ and x two-component Pauli spinors creating a heavy quark and a heavy 
antiquark, respectively, cr' is a Pauli matrix, and D is the spatial part of the gauge covariant 
derivative. Note that only a single polarization state is projected out in the expressions in 
Eqs. (12c) and (12d). In Eq. (12d), the following notation for the symmetric traceless 
component of a tensor is used: A^"^^^ = ^{A^^ + A^^) — ^5'^^ A'^''. The factors of 2E appear 
on the right sides of Eq. (12), because the state | [QQ]i{'^^~^^Lj)) on the left sides is normalized 
relativistically. In Eq. (8), we have taken the derivative with respect to q. Therefore, when 
we read the corresponding short-distance coefficients, we must cancel the overall factor |q| 
in Eqs. (12b), (12c), and (12d). Replacement of the gauge covariant derivative D by an 
ordinary derivative V brings in corrections of relative order v'^ in the Coulomb gauge [4] in 
which the NRQCD operators are defined. 

By making use of Eqs. (7), (9), and (12), we obtain the hadronic tensors: 



This step is equivalent to replacing the perturbative QQ state \[QQ]i{^^^^Lj)) in Eq. (11) 
by the meson state ^j2mH\H{^^^^Lj)) , where mn is the mass of H. At leading order in v, 
niH — 2£'|q=o = 2m. The meson state |if(^*+^Lj)), that is included in the vacuum-saturated 
analog of the NRQCD decay matrix elements for L = 5" or P, has the nonrelativistic 

normalization: 



|(//(%)|V'^x|0)r = (Oi)5, (14a) 

(i/(3Po)|^^(-|B-(T)x|0)|' = (Oi)p, (14b) 

I \{H('PiM^ (-|B XfT-en) x|0)f = (Oi)p, (14c) 

I J2{HeP2M^ (-|BV)6|) x|0)f = (Oi)p. (14d) 



The expressions in Eq. (14) have errors of relative order v^, which break the heavy-quark 
spin symmetry. 

Substituting Eq. (11) into Eq. (13), we obtain the hadronic tensor A^^ in Eq. (1): 



■^H(-^Pi) 



4e^e 



2^2 



m 



y3(g2 _ p2)2 
V2m(g2 _ p2)2 



2 



mi 



(15a) 
(15b) 



— I 



.^e'elPh*H,^{XH) 



(15c) 



X 



(Q2 - P^)^ 



g'^'Q'^Q'^ + P'^Q'^g"'^ - P'^Q'^g"'^ + -{Q^ - P^)g'^g'''^ 



(15d) 



The hadronic tensor (15) can be compared with previous results for the color-octet processes 
J + g* ^ [QQ]8{^^~^^Lj) [19]. We find that the results in Eq. (15) are consistent with those 
given in Eqs. (3.23), (3.25), (3.26), and (3.27) of Ref. [19] up to color factors. We have also 
checked that Eq. (15) reproduces the decay widths for rjQ, xqo, and xq2 into two photons [4] 
in the hmit — > 0. 

D. Cross section 



Substituting the hadronic tensor A'^ in Eq. (15) and polarization states for the H and 7 
into Eq. (1), we obtain the helicity amplitude. By squaring the helicity amplitude, averaging 
over the spins of the e'^ and e~ , dividing by the incident flux 2s, and multiplying by the 
two-body phase space, we obtain the differential cross section for each helicity state in terms 

of X = cos^^, where 9 is the scattering angle of the photon at the CM frame. The resultant 
differential cross section for the e'^e~ — > H{^So) + 7 process is 



da , 



e^e 



HCSo)+7{±l)] 



2 4 3 

TT Cqa r 



1 - 



m 



H 



where — 1 < x < 1 and the scaling parameter r is defined by 

4m^ 
r = . 



(l + x^) 



(16) 



(17) 



The differential cross section of the e+e~ — > Hi^Pj) + 7 process for each hehcity state is 
given by 



da . 



where nonvanishing entries of Fj{Xh, A) are 



(18) 



Fo(0,±l) = ^i-^(l + x2), (19a) 

Fi(0,±l) = l+x^ (19b) 

Fi(±l,±l) = 2r(l -x^), (19c) 

F2(0,±l) = + x^), (19d) 

F2(±l,±l) = 2r(l-x'), (19e) 

F2(±2,±l) = 2r'(l + a;2). (19f) 

In computing the phase space in Eqs. (16) and (18), we take the physical mass ttih instead 
of the invariant mass 2E of the QQ pair. This accounts for the factor 1 —m^/s in Eqs. (16) 
and (18). 

When the CM energy ^/s is much greater than the heavy-quark mass m, the asymptotic 
behavior of the cross section for e'^e" — > if + 7 satisfies the hehcity selection rules of 
perturbative QCD [20, 21]. In this hmit, the scaled cross section R{H + 7) in units of the 
cross section for approaches the asymptotic form [5] 

R[H{Xh) + 7] ~ o;t;='+'"(^V^+I^«I , (20) 

where n{S) — and n{P) — 1. The factor ■i;3+2n(L) -^^ ^20"^ arises from the wave 
function of the H. One power of r in Eq. (20) stems from the large momentum transfer 
which is required for the heavy-quark pair to form the heavy quarkonium with small relative 
momentum and the other powers arise from the helicity selection rules [5]. In the limit 
r — > 0, the differential cross sections in Eqs. (16) and (18) satisfy the asymptotic form (20). 
For the exclusive charmonium production associated with a photon at the B factories, the 
ratio r fa 0.07 is small for rric — 1.4 GeV. However, for the case of a bottomonium, r fs 0.76 
is of order 1 for irif, — 4.6 GeV. Therefore, the asymptotic behavior (20) is approximately 
satisfied only for the exclusive charmonium plus 7 production at the B factories. 
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Integrating the differential cross section given in Eq. (16) over x and summing over the 
photon hehcities, we find the total cross section for the process e+e"~ — > H{^Sq) + 7: 

,[,.e-^i/(X)+7l = ^(l-^)<^^. (21) 

6s \ s J m"* 

Integrating Eq. (18) over x and summing over the helicities of the photon and quarkonium, 
we obtain the total cross sections for e^e~ — > Hi^Pj) + 7: 

9s(l-r)^ V s / ' 

where we have used the spin-1 and spin-2 polarization tensors for ^Pi and states summed 
over polarization states: 

Xh 

^HMe*H'^{XH) = lil^^rl^ + I'^^F'^) - i/z^-r^. (23b) 

Xh 2 3 

III. PARAMETERS FOR THE NUMERICAL ANALYSIS 

In order to predict the cross sections for e^e~ t/q + j and xqj + 7 at the B factories 
by using Eqs. (21) and (22), we need to determine the numerical values for the NRQCD 
matrix elements and input parameters such as m, mn, ck, and ag. In this section, we 

list the numerical values for those quantities that we use in this work. 



A. Input parameters 

The short-distance coefficients for the cross sections in Eqs. (21) and (22) depend on the 
heavy-quark mass m. We take the one-loop pole mass for the heavy-quark mass m: rric = 
1.4±0.2 GeV for the charm quark and — 4.6±0.1 GeV for the bottom quark, respectively. 
As we mentioned earlier, we use the physical quarkonium mass for niH in the phase-space 
factor 1 - mjj/s in Eqs. (21) and (22), where m^^ = 2.9804 GeV, m^^(25) = 3.638 GeV, 
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m^^o = 3.41475 GeV, m^^, = 3.51066 GeV, m^^^ = 3.55620 GeV, m^^ = 9.585944 GeV, 
m^^^ = 9.89278 GeV, and m^^^ = 9.991221 GeV [22]. For the rjb mass, we take m^, = 9.3889 
GeV [23]. 2 

The cross sections in Eqs. (21) and (22) are proportional to a^. For the QED couphngs 
for the virtual photon, we use the running coupling constant [24] = ^/s) = 1/131 at 
— 10.58 GeV. The QED couphng for the real photon is chosen to be q;(/x = ^) w 1/132, 
where z — [{s — 2m^)/2]^/^ is the invariant mass of the virtual heavy quark: z 7.3 GeV 
and 5.9 GeV for the charmonium and the bottomonium, respectively. 

B. Long-distance matrix elements 

1. NRQCD matrix elements {Oi)s 

The color-singlet NRQCD matrix elements for the S'-wavc spin-singlct charmonia, rjc and 
r]c{2S), can be determined from the two-photon decay rates. However, up to errors of relative 
order v"^ that arise from the heavy-quark-spin-symmetry breaking, one can determine the 
values also by making use of the leptonic decay rates of the spin-triplet counterparts, which 
are far more accurately measured. Among various theoretical attempts to determine the 
matrix elements for the -S'-wave spin-singlet charmonia, we employ a recently developed 
method given in Refs. [14, 25, 26]. This method resums a class of relativistic corrections to 
the color-singlet contributions to the electromagnetic decay rates of both spin-singlet and 
triplet 5"- wave charmonia to all orders in v. The color-singlet NRQCD matrix element {Oi)s 
for the 77c is taken to be OASltoHl GeV^ [26]. 

For the radially excited S-wme spin-singlet state r}c{2S), the NRQCD matrix element 
determined from the measured decay width for rjci^S) 77 differs significantly by about 
4.60" from that for ip{2S) determined from ip{2S) — > e'^e'^ . The difference is greater than 
the errors of neglecting spin-symmetry breaking effects. We suspect that the main reason 
for the discrepancy is originated from the assumption [27]: 

Br[77e(25) ^ K^K^TT^] = Br[rj, ^ K^s^^tt^]. (24) 

^ Very recently, the BABAR Collaboration observed the rjb resonance in the photon energy spectrum in 
radiative T(3S') decay [23]. 
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This assumption was imposed on the indirect determination of the partial width for 
r)c{2S) 77 from the measurement for 7*7* r)c{2S) KgR^n^ [27]. In this work, 
instead of using this indirect determination of F [77c (25') 77] in Ref. [27], we assume that 
the NRQCD matrix element determined from r[ric{2S) — > 77] must be the same as that de- 
termined from r[iJj{2S) — ^ e~^e~] within errors of relative order f^, which breaks the heavy- 
quark spin symmetry. Our assumption results in the prediction for the two-photon width of 
ric{2S): r[77c(25') — > 77] = 2.92 keV. Direct measurement of the two-photon width of rjc{2S) 
may provide one with a stringent test of our assumption. We will present a detailed analysis 
regarding this point in a separate publication. Basic strategy for determining the matrix ele- 
ment is similar to that in Ref. [26]. Disregarding the assumption (24) and making use of the 
matrix element obtained from r[ilj{2S) e^e~] based on the heavy-quark spin symmetry 
up to corrections of relative order v^, we find that {Oi)s = 0.274+|j:236 GeV^ for S = r]c{2S). 

The NRQCD matrix element {Oi)s for the 77b cannot be determined from 77^ — > 77 because 
the measurement is not available, yet. Instead, we quote the value for the matrix element 
(Ci)t(is) given in Eq. (23a) of Ref. [28], which is determined from T(15') e+e~. We 
assume that the heavy-quark spin symmetry is relatively well satisfied because ~ 0.1 is 
small for the bottomonium. The numerical value is {Oi)s = 3.069lo;^gQ GeV^ for S = r]h [28]. 



2. NRQCD matrix elements {Oi)p 

The generalization of the resummation of relativistic corrections [14] to the P-wave 
quarkonium has not been developed because of the complication due to the color-octet 
contributions. Therefore, for the spin-triplet P-wave charmonium XcJ, we determine the 

NRQCD matrix element {Oi)p by comparing the f-leading NRQCD factorization formulas 
for the photonic decay widths [4] 



r(Xco 77) = 67re^a m 



1 + 



(Stt^ - 28)q;, 



187r 



2 



{O 



IP 



(25a) 



rfe-TT)^^^l-^j ^ (25b) 

with corresponding measured values. The decay of Xci into 77 is absent due to Yang's 
theorem [29]. 

Substituting the numerical values Br(xco 77) — (2.76±0.33) x 10~^ [22] , the total decay 
width r^^o = 10.4 ± 0.7 MeV [22], ^^(/x = m^^j/2) = 0.32, and a{n = m^^^/2) = 1/133 [24] 
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into Eq. (25a), we obtain (Oi)p = {Oi)]^^^ = 0.051 ± 0.010 GeV^ for P = Xco, where the 
superscript 77 indicates that the value is fit to the two-photon width. Taking into account 
Br(xc2 ^ 77) = (2.58 ± 0.19) x 10"^ [22] and T^^^ = 2.05 ± 0.12 MeV [22], we obtain the 
matrix element (Ci)p = {^1)12,2 ~ 0-068 ± 0.009 GcV^ for P = Xc2- Here, the uncertainties 
from the heavy-quark mass are not included in both (C^i)^2o (^i)xl2' t^]' value 

(^i)xc2 — 0.053 ± 0.009 GeV^ was obtained by making use of a similar method. The main 
reason for the difference is the change in the measured branching fraction of the Xc2 decay into 
77 from (2.19 ± 0.33) x 10"^ [30] to (2.58 ± 0.19) x 10"^ [22]. We take the weighted average 
over the spin states of Xco and Xc2 to determine (Oi)p = 0.060^0 029 GeV^ for P — XcJ, 
where the heavy-quark spin symmetry is imposed and the uncertainty from the heavy-quark 
mass is considered together with those from the experimental values. 

In the cases of the spin-triplet P-wave bottomonium states X5j(lP), the measured 
values for the two-photon decay widths are not available. Therefore, we use the first 
derivative of the wave function at the origin that has been determined by using the 
Buchmiiller-Tye potential [31]. The resultant value for the NRQCD matrix element is 
(Ci)p = 2.03 GeV^ for P = XbJ [32]. 

3. Summary of NRQCD matrix elements 

In Table I, we tabulate the color-singlet NRQCD matrix elements for rjc, r]c{2S), rjb, 
XcJ-i and XhJ- As has been described in the text, there are input parameters such as m 
that have been used to determine the NRQCD matrix elements. When we compute the 
short-distance coefficients for the cross section for e^e~ H + ^, such input parameters 
must be consistent with the values that have been used to determine the NRQCD matrix 
elements [26]. In order to take into account such correlations, we present sources of the 
uncertainties for each NRQCD matrix element in Table I. The first and the second lines 
after headings contain uncertainties arising from the heavy-quark mass m and from sources 
described in each reference except for m, respectively. The sum of the two uncertainties is 
obtained by adding them in quadrature as shown in the last row of Table I. 
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TABLE I: NRQCD matrix elements {Oi)h for H{^So) in units of GeV^ and for H{^P,j) in units of 
GeV^. The first and the second Hnes after headings contain uncertainties arising from the heavy- 
quark mass m and from other sources described in the text. Uncertainties in the last line are 
obtained by adding the two uncertainties in quadrature. 



Sources of errors \ H 


Vc [26] 


Vc{'2S) 


% [32] 


XcJ 


Xbj [31, 32] 




(GeV3) 


(GeV3) 


(GeV3) 


(GeV^) 


(GeV^) 


Am 


r. 407+0.033 

^■^"5' -0.025 


'-'■^'^-0.010 


n«Q+0-000 

O.UU»_Q QQQ 


U.UDU_Q 


2.03 


others 


Q 407+0. 106 


Q 274+0-040 


O.UU»_Q^gO 


0.060 ± 0.007 


2.03 


total 


Q 407+0.111 


Q 274+0-042 

'-'■^'^-0.036 


O.UU»_QjgQ 


060+° °'^-^ 

U.UDU_Q 


2.03 



IV. PREDICTIONS FOR THE B FACTORIES 



In this section, we provide the predictions for the cross sections for e^e^ ^ + 7 at 
the B factories. We also discuss the uncertainties and phenomenological imphcations of the 
results. 

Our predictions for the total cross section a for e^e" — > if + 7 at the CM energy 
y/s — 10.58 GeV are given in the first column of Table II. The values have been obtained by 
substituting the input parameters listed in Sec. Ill into the NRQCD factorization formulas 
in Eqs. (21) and (22). In the second column, we provide the cross section a^nt in which we 
have imposed the cut | cos^j < 0.8 for the scattering angle of the photon at the CM frame. 
(Tcut is about 70%-80% of a in any case. In the last column in Table II, we list the energy 

— [s — rriff) I i^\fs) of the photon emitted in e^e" — > if + 7 at the CM frame, which may 
provide an efficient trigger for the signals. 

A. Uncertainties 

The error bars of the cross sections a and o"cut listed in Table II contain the following 
uncertainties. A dominant source of the uncertainties is the heavy-quark mass m, which 
affects both the NRQCD matrix elements in Table I and the short-distance coefficients 
appearing in the NRQCD factorization formulas in Eqs. (21) and (22). To avoid double 
counting of those errors, we have considered the correlation of the uncertainties from m. 
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TABLE II: The total cross section a in units of fb and the photon energy at the CM frame in 
units of GcV in e+e~ H + -f depending on charge-conjugation parity C = +1 quarkonium H. 
The subscript "cut" represents the cut | cos^| < 0.8 for the scattering angle 9 of the photon at the 
CM frame. 



H 




<7"cut (fb) 


(GeV) 


Vc 




59.7tlli 


4.87 


Vc{2S) 


49.2l?;4 




4.66 


rib 


o c+0.2 
^■"^-0.2 


i.8i°:l 


1.12 


XcO 


1 q+0.2 

^■''-0.2 


1.01°:} 


4.74 


Xcl 


13.7l|t 


10.211^ 


4.71 


Xc2 






4.69 


XbO 


U.D_o.2 


0.41°:? 


0.95 


Xbl 




9 q+0.6 


0.66 


Xb2 




2 4+0.8 

^•^-0.5 


0.57 



Other errors in {Oi)s include a theoretical uncertainty of relative order from the leading- 
potential approximation, the uncertainty of the string tension, which is a parameter of 
the Cornell potential model, the uncertainty reflecting the uncalculated next-to-next-to- 
leading-order corrections in as to the electromagnetic decay width of the quarkonium, and 
the uncertainty of the measured decay width. For more details about the error analysis, we 
refer the reader to Refs. [26]. The error bars of the NRQCD matrix elements {Oi)p for the 
Xcj in Table I include the uncertainties from the heavy-quark mass m and the measured 
widths for the two-photon decays of the Xco and Xc2- The NRQCD matrix elements {Oi)p 
for the XbJ a-re quoted from Refs. [31, 32], where uncertainties are not included [32]. As a 
result, except for the uncertainties from the heavy-quark mass m appearing in the short- 
distance coefficients, no additional uncertainties are included in the error bars of the cross 
section for H — XbJ- The error bars of a and Ucut in Table II have been obtained by adding 
all the uncertainties listed above in quadrature. 

Our predictions in Table II may suffer from QCD and relativistic corrections that have 
not been calculated, yet. We guess that the corrections may not be as large as those observed 
in the two-charmonium production from e'^e" annihilation [12, 15, 16, 33], where all such 
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corrections are aligned to make a significant enhancement of the prediction for the cross 
section. 

B. Exclusive r]c{lS, 2S) + 7 production 

As we have shown in Table 11, the cross sections for the S'-wave spin-singlet charmonia 
H = rjc and rjci'^S) produced exclusively with a photon are about 80 fb and about 50 fb, 
respectively. These values are significantly greater than those for J/i/j + rjc and J/jlj + rjc{2S) 
measured at the B factories. We have predicted that r[77c(25') — > 77] = 2.92 keV by replacing 
the assumption (24) that was imposed on the indirect determination of r[77c(2S') — > 77] 
in Refs. [34, 35] with an alternative assumption that the approximate heavy-quark spin 
symmetry still holds between the NRQCD matrix elements for r]c{2S) and iIj{2S). The 
measurement of the cross section for e~^e~ — > ric{2S) -|- 7 is particularly interesting because 
the measurement may provide an independent test of our argument presented in Sec. IIIB 1. 

C. Exclusive 775-1-7 production 

Our prediction of the cross section for H — rj^ a,t the B factories is about 2.5 fb. By 
taking into account the integrated luminosities available at the present B factories, we 
expect that over three thousands of e"^e~ — > + 7 events can be produced. This number 
exceeds the estimated number of events at LEP II through 7*7* rjb by about an order of 
magnitude [36-38]. Therefore, the exclusive process e'^e~ — > 776 + 7 could be an alternative 
probe to the rjf, meson as long as the background can be removed significantly. 

D. Exclusive xqj + 7 production 

As shown in Table II, the cross section for H = xqj ranges from about 1 to 5 fb except 
for the case of Xci, which has a particularly large value of about 14 fb. The reason is that 
the overall coefficient for the longitudinal xqi is greater than those for the other two states 
as shown in Eq. (19).^ This contrasts to the P-wave decay into 77, where XQi ~^ 77 is 
forbidden by Yang's theorem [29]. A similar feature can be found in the inclusive charm 

^ According to the hadron helicity selection rules, the longitudinal cross section for XQJ dominates as r — > 0. 



17 



production in XbJ decay [32]: the decay rate of the Xbi into g*g is comparable to those of 
the Xbo and the Xb2, especially in the limit of the large invariant mass of the g*. Unlike the 
P-wave spin-triplet charmonium, the cross sections for Xbi a-nd Xb2 are comparable and that 
for Xbo is much suppressed. 



E. Backgrounds 

The exclusive process e'^e~ — > i?-|-7 for C = +1 heavy quarkonium H suffers from large 
contamination from the background e'^e" — > X + 7. In this section, we provide a rough 
estimate on the background. 

We assume that the dominant subprocesses of the background e'^e' X + j are e'^e" 
qq + 7, where q = u, d, s, and c. These processes have collinear divergences due to the 
photon radiation from the initial leptons and final massless quarks. In order to avoid such 
complications, we take the following cuts: |cos^| < 0.8, cos9qj < 0.8, and cos^^^ < 0.8, 
where 9 is the scattering angle of the photon and is the relative angle between the photon 
and a quark or an antiquark i in the CM frame. We also restrict the energy range of the 
photon by requiring the recoil mass mx to satisfy \mx — ^//| < 50 MeV. 

By imposing the constraints listed above, we carry out the numerical calculation for the 
background cross section aBciH) for each C = +1 heavy quarkonium H by making use of 
CompHEP [39]. The results are (Jbg (charmonium) « 0.5 pb and crBG(bottomonium) Ri 2.4- 
4.5 pb. The background cross sections aBciH) are much greater than the signal cross 
sections dcut in Table II. However, the high integrated luminosities C of order 1 ab"^ of 
the present B factories might allow one to measure the exclusive C = +1 quarkonium 
production associated with a photon. We can compute the signal significance S{H) for each 
heavy quarkonium, where S{H) is defined by 

s(H) = "7'^+;;'f . (26) 

We find that 5*- wave charmonium states have particularly large values: S{r]c) ~ 81 and 
S[ric{2S)] ~ 47. In the cases of P-wave charmonium states, <S'(xco) ~ 1) S{xci) ~ 13, and 
S{Xc2) ~ 5. And S{H) is less than 1 for rjb and XbJ states. 

Our background study given above is only a rough estimate. In the cases of bottomo- 
nium states, there must be obstacles in detecting photons with energies about 1 GeV or 
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less because of the large contamination from the decays of hght hadrons [40]. However, 
the exclusive charmonium production associated with a photon may have significantly less 
additional backgrounds unlike the bottomonium states because the radiating photon is suffi- 
ciently hard. We also note that there may be another background arising from the fecddown 
from the e"'"e^ H(^Si) + 7 process for H{^Si) = J/ip, ip{2S), or T{nS), whose cross sec- 
tions are of order pb. Especially, the measurement of the XQj-l-7 production may be affected 
by such feeddown. Reduction of such backgrounds may be subject to detector resolution. 

V. SUMMARY 

We have calculated the cross section for the C = +1 heavy quarkonium H produced 
exclusively with a photon from e+e^ annihilation into a virtual photon at leading order in 
as and v. The NRQCD factorization formulas for the differential distributions with respect 
to the scattering angle of the photon and the total cross sections for H — rjQ and xqj have 
been obtained in Eqs. (16), (18), (21), and (22), respectively. Our predictions for the cross 
sections at ^/s — 10.58 GeV are tabulated in Table II. The predicted cross sections for the 
5"- wave spin-singlet states are about 80 fb, 50 fb, and 3 fb for rjc, ric{2S), and rjh, respectively. 
In the cases of the P-wavc spin-triplet charmonia, the cross section for Xci is about 14 fb; 
that is significantly greater than those for Xco and Xc2- The cross sections for Xbi and Xb2 
are about 3 fb and that for x^o is suppressed. Under the cut | cos^| < 0.8 for the scattering 
angle 9 of the photon, the cross sections decrease by about 20%-30 %. 

The exclusive process e+e~ — > if -|- 7 may suffer from a large background from e'^e~ — > 
X + J. A rough estimate on the background cross section has been made by considering 
only the e^e~ ^ qq + '~f contributions at leading order in a^. In computing the background 
cross section, we have required \mx — ^h\ < 50 MeV for the recoil mass mx and imposed 
the angular cuts |cos^| < 0.8, cos^^-y < 0.8, and cos^^-y < 0.8. The resultant background 
cross sections are about 0.5 pb for the charmonia and about several pb for the bottomonia. 
The signal significances for the rjc, r]c{2S), Xci, and Xc2 cases are greater than 6. However, 
those for the rjb, Xco, and XbJ are about 1 at the present B factories. 

Very recently, the BABAR Collaboration observed the rih resonance in the photon energy 
spectrum in radiative T(3S') decay [23]. The exclusive process e~^e~ ^ f]b + 'y that we have 
considered could be an alternative probe to the rji, meson as long as the backgrounds are 
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significantly removed. 

As we have remarked earlier, the authors of Refs. [34] and [35] made an assumption 
Br[r/c(25) ^ KIK^-k"^] = Br[r/c -> K'^gK^n^] in determining the T[r]^{2S) 77] = 1.3 ± 
0.6 keV [34] and 0.59 ± 0.19 keV [35], respectively. This assumption seems to be the origin 
of inconsistency of the NRQCD matrix element for r]c{2S), which was fit to r[r]c{2S) — >■ 
77] [34, 35], in comparison with that for ip{2S). As an alternative choice, we assume that 
the approximate heavy-quark spin symmetry holds between the NRQCD matrix elements 
for the r]c{2S) and ^(2-5) to predict r[r]c{2S) 77] = 2.92 keV. The measurement of the 
cross sections for e'^e" Vc+J and e'^e~ —>■ r]c{2S)+j may provide a stringent independent 
constraint to test our argument. 
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